LOFAR sparse image reconstruction

The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SK

Segmentation of Loops from Coronal EUV Images

We present a procedure which extracts bright loop features from solar EUV images. In terms of image intensities, these features are elongated ridge-like intensity maxima. To discriminate the maxima, we need information about the spatial derivatives of the image intensity. Commonly, the derivative estimates are strongly affected by image noise. We therefore use a regularized estimation of the derivative which is then used to interpolate a discrete vector field of ridge points ``ridgels'' which are positioned on the ridge center and have the intrinsic orientation of the local ridge direction. A scheme is proposed to connect ridgels to smooth, spline-represented curves which fit the observed loops. Finally, a half-automated user interface allows one to merge or split, eliminate or select loop fits obtained form the above procedure. In this paper we apply our tool to one of the first EUV images observed by the SECCHI instrument onboard the recently launched STEREO spacecraft. We compare the extracted loops with projected field lines computed from almost-simultaneously-taken magnetograms measured by the SOHO/MDI Doppler imager. The field lines were calculated using a linear force-free field model. This comparison allows one to verify faint and spurious loop connections produced by our segmentation tool and it also helps to prove the quality of the magnetic-field model where well-identified loop structures comply with field-line projections. We also discuss further potential applications of our tool such as loop oscillations and stereoscopy.Comment: 13 pages, 9 figures, Solar Physics, online firs

Multiresolution analysis of active region magnetic structure and its correlation with the Mt. Wilson classification and flaring activity

Two different multi-resolution analyses are used to decompose the structure of active region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average and standard deviation of the magnetic flux gradient for alpha, beta, beta-gamma and beta-gamma-delta active regions increase in the order listed, and that the order is maintained over all length-scales. This study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length-scales in the active region, and not just those length-scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length-scales. It is also shown that the average gradient content of active regions that have large flares (GOES class 'M' and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length-scales.Comment: Accepted for publication in Solar Physic

Templates for Convex Cone Problems with Applications to Sparse Signal Recovery

This paper develops a general framework for solving a variety of convex cone problems that frequently arise in signal processing, machine learning, statistics, and other fields. The approach works as follows: first, determine a conic formulation of the problem; second, determine its dual; third, apply smoothing; and fourth, solve using an optimal first-order method. A merit of this approach is its flexibility: for example, all compressed sensing problems can be solved via this approach. These include models with objective functionals such as the total-variation norm, ||Wx||_1 where W is arbitrary, or a combination thereof. In addition, the paper also introduces a number of technical contributions such as a novel continuation scheme, a novel approach for controlling the step size, and some new results showing that the smooth and unsmoothed problems are sometimes formally equivalent. Combined with our framework, these lead to novel, stable and computationally efficient algorithms. For instance, our general implementation is competitive with state-of-the-art methods for solving intensively studied problems such as the LASSO. Further, numerical experiments show that one can solve the Dantzig selector problem, for which no efficient large-scale solvers exist, in a few hundred iterations. Finally, the paper is accompanied with a software release. This software is not a single, monolithic solver; rather, it is a suite of programs and routines designed to serve as building blocks for constructing complete algorithms.Comment: The TFOCS software is available at http://tfocs.stanford.edu This version has updated reference

Power, norms and institutional change in the European Union: the protection of the free movement of goods

How do institutions of the European Union change? Using an institutionalist approach, this article highlights the interplay between power, cognitive limits, and the normative order that underpins institutional settings and assesses their impact upon the process of institutional change. Empirical evidence from recent attempts to reinforce the protection of the free movement of goods in the EU suggests that, under conditions of uncertainty, actors with ambiguous preferences assess attempts at institutional change on the basis of the historically defined normative order which holds a given institutional structure together. Hence, path dependent and incremental change occurs even when more ambitious and functionally superior proposals are on offer

Population overlap and habitat segregation in wintering Black-tailed Godwits Limosa limosa

Distinct breeding populations of migratory species may overlap both spatially and temporally, but differ in patterns of habitat use. This has important implications for population monitoring and conservation. To quantify the extent to which two distinct breeding populations of a migratory shorebird, the Black-tailed Godwit Limosa limosa, overlap spatially, temporally and in their use of different habitats during winter. We use mid-winter counts between 1990 and 2001 to identify the most important sites in Iberia for Black-tailed Godwits. Monthly surveys of estuarine mudflats and rice-fields at one major site, the Tejo estuary in Portugal in 2005-2007, together with detailed tracking of colour-ringed individuals, are used to explore patterns of habitat use and segregation of the Icelandic subspecies L. l. islandica and the nominate continental subspecies L. l. limosa. In the period 1990-2001, over 66 000 Black-tailed Godwits were counted on average in Iberia during mid-winter (January), of which 80% occurred at just four sites: Tejo and Sado lower basins in Portugal, and Coto Dontildeana and Ebro Delta in Spain. Icelandic Black-tailed Godwits are present throughout the winter and forage primarily in estuarine habitats. Continental Black-tailed Godwits are present from December to March and primarily use rice-fields. Iberia supports about 30% of the Icelandic population in winter and most of the continental population during spring passage. While the Icelandic population is currently increasing, the continental population is declining rapidly. Although the estuarine habitats used by Icelandic godwits are largely protected as Natura 2000 sites, the habitat segregation means that conservation actions for the decreasing numbers of continental godwits should focus on protection of rice-fields and re-establishment of freshwater wetlands

LOFAR sparse image reconstruction

International audienceContext. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods. Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the " compressed sensing " (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework. Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data. Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2−3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A-and W-projections) required for current and future instruments such as LOFAR and SKA

Euclid: On the reduced shear approximation and magnification bias for Stage IV cosmic shear experiments

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculate the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, n_s, sigma_8, Omega_DE, w_0, and w_a of -0.51 sigma, -0.36 sigma, 0.37 sigma, 1.36 sigma, -0.66 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant